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We establish a method to trace the Lagrangian evolution of extended objects consisting of a mul- 
ticomponent scalar field in terms of a numerical calculation of field equations in three dimensional 
Eulerian meshes. We apply our method to the cosmological evolution of global strings and evaluate 
the energy density, peculiar velocity, Lorentz factor, formation rate of loops, and emission rate of 
Nambu-Goldstone (NG) bosons. We confirm the scaling behavior with a number of long strings per 
horizon volume smaller than the case of local strings by a factor of ~ 10. The strategy and the 
method established here are applicable to a variety of fields in physics. 



PACS number(s): 98.80.Cq 

Extended objects consisting of scalar fields such as 
topological defects or solitons play important roles in 
many fields of physical science ranging from condensed 
matter physics to cosmology. Much work has been done 
on the cosmological formation of topological defects by 
thermal or nonthermal |^ phase transitions, their sub- 
sequent evolution ||^ , and cosmological implications of Q- 
balls On the other hand, it has been pointed out that 
cosmological defect formation may be realizable in a lab- 
oratory ||5|] . The phase transition of liquid helium can also 
be described in terms of Ginzburg-Landau theory as in 
the case of cosmological phase transitions. The order pa- 
rameter, which is an expectation value of a scalar field in 
cosmological phase transition, is just replaced by a Bose 
condensate wave function. Following Zurek's suggestion, 
the formation of topological defects has been studied in 
■^He Q and ''He j^]. In fact, vortex formation is observed 
in '^He. Thus, it is very important to investigate the 
formation and subsequent evolution of extended objects, 
such as topological defects, produced both in cosmology 
and laboratories. 

Among various types of topological defects predicted 
in high energy physics, strings hold a unique position, 
which, unlike magnetic monopoles and domain walls, do 
not overdose the universe because they settle down to 
the scaling solution The key mechanism to achieve 

it is intercommutation of infinite strings to lose their en- 
ergy by producing closed loops, which decay-through ra- 
diating gravitational waves for local stringsD or Nambu- 
Goldstone (NG) bosons for global strings. 

While many grand unified theories may contain lo- 
cal strings, there are no particle physics motivations for 
their existence. Nonetheless they were extensively stud- 
ied because it was expected that they might be the origin 
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of density fluctuations which seeded cosmic large-scale 
structures. In particular, time evolution of local cosmic 
strings have been investigated by a number of authors 
[|l^,[l]jl2[|l|] who used the Nambu-Goto action @ that 
applies for infinitely thin strings. These numerical anal- 
yses confirmed the scaling behavior and the scaling pa- 
rameter, ^, which is defined as ^ = pst^ / ^J, with ps energy 
density of strings and p string tension per unit length, 
converged to ^ ~ 10 in the radiation dominated universe 



*The possibility has also been pointed out that infinite 
strings lose their filiergy by directly emitting massive particles 
for local strings |9|. 



13 1 . Recent observations of anisotropics of cos- 
mic microwave background radiation, however, disfavor 
cosmic-string scenario of structure formation |l5| , and 
motivations to study local strings have somewhat dimin- 
ished. 

Global strings, on the contrary, are well motivated 
in particle physics, whose existence is predicted by the 
Peccei-Quinn solution to the strong CP problem [ p^ . 
They are produced when the Peccei-Quinn U(l) symme- 
try is broke n ||l6| , p7[ |. They radiate axions as associated 
NG bosons |18| which are one of the best candidates of 
cold dark matter. Furthermore, they are also realized 
in condensed matter physics. Vortices of ■^He and ^He 
studied in a laboratory are analogous to global strings 

In contrast to local strings, energy density of global 
strings is dominated by the gradient energy of the NG 
scalar field rather than the potential energy of the core, 
because there are no gauge fields to cancel the former. As 
a result they have a long-range interaction and their dy- 
namics cannot be analyzed with the Nambu-Goto action. 
It should be described with the Kalb-Ramond action in- 
stead , which incorporates NG bosons and their cou- 
pling with the string core in addition to the Nambu-Goto 
action. In this action, however, the line density of the 
string is not well defined which has a logarithmic depen- 
dence on the radius. It is also very difficult to analyze it 
numerically. Hence the dynamics of strings governed by 
the Kalb-Ramond action has not been clarified in cos- 
mological context yet even though evolution of axionic 
strings has very important implications in cosmology and 



1 



particle physics. Instead, the results of numerical simu- 
lations of local strings have been used to discuss cosmol- 
ogy of global strings as well, which may not be justified 
because evolution of global defects can be significantly 
different from that of local counterparts as seen in the 
case of local and global ||2^,|2^ monopoles. 

In order to surmount such a situation, we have inves- 
tigated evolution of global strings by solving equation of 
motion of the complex scalar field numerically in a three 
dimensional lattice in the expanding universe instead of 
dealing with the Kalb-Ramond action in collaboration 
with Kawasaki |2jj2^. In these fully Eulerian simula- 
tions, it is difficult to identify strings and follow their 
dynamics. In the previous publications, whether a lat- 
tice contains a string core was judged by comparing po- 
tential energy density at each vertex with a static cylin- 
drically symmetric solution of a global string. Counting 
number of lattices thus identified as containing a string, 
we estimated the scaling parameter ^ and found that 
f ~ 1, which is significantly smaller than the case of local 
strings. This result has been criticized by Martins and 
Shellard |^,^ who claim that we have obtained incor- 
rectly smaller value of the scaling parameter because the 
energy loss of long strings due to direct emission of NG 
bosons is much larger in simulation with small dynamic 
range than in the actual cosmological context, and that 
both types of long strings should behave quantitatively 
in the same way on cosmological scales. 

In order to characterize cosmological evolution of 
global strings completely, we need to establish the 
method to follow cosmological evolution of global strings 
from a Lagrangian point of view. The energy density of 
long strings, poo satisfies the following equation: 



dpo 



dt 



-2H{1 + {v^))poo - TloopPoo - TngPc 



(1) 



where the second and the third terms on the right-hand- 
side represent energy loss due to loop formation and di- 
rect emission of NG bosons or axions, respectively, and 
(w^) denotes average square velocity of string segments. 
In the scaling regime, the string network is characterized 
by a scale L, which is defined as L = y/VPoo ^-^^ grows 
with the horizon scale L (x t. Here /i is the effective line 
energy density of a global string. We can introduce a 
loop production coefficient c and an emission coefficient 
K by 



Poo ^ Poo 

i loopPoo = ^ NGPoo = 



(2) 



Then, the scaling parameter ^ is characterized by these 
coefficients as 



C + K 



(3) 



Indeed if k incorrectly turned out to be much larger than 
c, we would find a smaller value of ^ than it should really 



be. Hence it is essential to evaluate contribution of each 
term in (|^) with the help of numerical analysis. For this 
purpose we must know how each string segment moves 
and intercommutes with each other to trace formation of 
loops, which is impossible in conventional Eulerian sim- 
ulations. Notice that this problem is ubiquitous because 
in many fields of physics we must solve equation of mo- 
tion in Eulerian meshes, for example, in order to trace 
time evolution of solitons like Q-balls [Q. 

In the present Rapid Communication, we establish a 
method to follow Lagrangian evolution of extended ob- 
jects in the context of Eulerian calculation of scalar field 
equations on the lattice. As a specific example, we con- 
centrate on global strings here, but the strategy and the 
method established here directly apply to other situa- 
tions in physics with some proper modifications. 

We express the complex scalar field, which constitute 
global strings, in terms of two real scalar fields (j}a{x) 
(a = 1, 2) and define their Lagrangian density as 
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(4) 



where g^i, is the fiat Robertson- Walker metric. The 
finite-temperature effective potential V[(j)ai T] is taken as 

V[cl>a.T]^]{cj^^-a^f + ^T^\ cl^'^^l + ct^l (5) 
4 6 

which exhibits a second-order phase transition with the 
critical temperature — -s/Sct which produces global 
strings by breaking a global U(l) symmetry. 

Equations of motion for the scalar fields are given by 



4>a{x) + 3H(j)a{x) - 



Rity- 



dV 



(6) 



where a dot denotes time derivative and R{t) is the cos- 
mic scale factor. In the radiation dominated universe, 
the Hubble parameter H = R{t)/R{t) and cosmic time t 
are given by 



3M^i 30 



(7) 



respectively, where Mp\ — 1.2 x lO^^GeV is the Planck 
mass and (7* is the total number of degrees of freedom for 
the relativistic particles. For the ease of numerical calcu- 
lations we take a = 0.1{45/16n^g^y/'^Mpi and A = 0.08 
but the results are insensitive to these choices. We start 
numerical simulation at the temperature Ti = 2Tc corre- 
sponding to U = tc/4: and adopt as an initial condition 
the thermal equilibrium state with a mass squared equal 
to the inverse curvature of the potential at that time. 
We have simulated the system from ten different thermal 
initial conditions under the periodic boundary condition. 
The number of lattices is 128"^ with the lattice spacing 
6x = ^/3tti/8 and the time step St = O.Olti. Thus, the 
box size is nearly equal to horizon volume (H^^)^ and 
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the lattice spacing to a core size of a string 5 = 1/ (v 2Act) 
at the final time t/ = 200ti. 

The above setup of calculation is the same as our pre- 
vious fully Eulerian simulations |2j^^^. In order to re- 
produce Lagrangian evolution of string segments in the 
Eulerian meshes, we develop a different method to iden- 
tify string cores from scalar field configurations, because 
the previous method [p4 25 1 was inadequate for config- 



urations with large curvature such as small loops and 
it was impossible to find more correct position of string 
core in a box beyond the lattice spacing. Here we use a 
new two-step algorithm to identify the locus of a string 
more accurately. The first step is to find plaquettes 
which a string penetrates. This is done in terms of the 
Vachaspati-Vilenkin algorithm ||2^ by monitoring phase 
rotation around each square. Then using the value of (pa 
at each vertex of a plaquette penetrated by a string, we 
linearly interpolate (j)a{x) to calculate the position where 
both (piix) and 4'2{x) vanish. We can thus find where a 
string penetrate in each plaquette and find more accurate 
trajectory of a string by connecting these points. 

The next and the more important task to obtain La- 
grangian evolution of string network, that has been miss- 
ing in the previous Eulerian analyses, is to find velocity of 
each string segment. Since motion tangential to a string 
is a gauge mode, we should evaluate velocity normal to 
it. Suppose that a string core exists at a point Xq at 
time toj namely, (/)i(a;o,to) — 02(3^0,^0) = 0. In order to 
estimate where this string segment moves At later, we 
expand scalar fields ipaiXjtQ + At) around ipaixojto) up 
to the first order, 

(j)a{x, to + At) = (paixo, to) + V 4>a{xo-, h) ' {x - Xo) 

+Mxa,to)At (a = 1,2). (8) 

Again the loci of string core a,t t = to + At are given by 
the line where both scalar fields (t)a{x,to + At) vanish. 
That is, under the first-order approximation (H) it lies on 
the intersection of two planes given by 

V0q(cco, to) • {x - Xo) + (}>a{xo, to) At = 0, (9) 

with a = 1,2. Now suppose that the line normal to 
the string segment at {x, to) intersect with the trajectory 
of string at t = to -|- At at a; = a;;(a;o, to. At). Then, 
making use of the fact that xi satisfies Eq. (^), we can 
easily calculate the velocity of this string segment, and 
its magnitude is given by 



a;;(£Co,to, At) - xq 



v{xo,to) = lim 

At >0 At 



|V(/)i X V(/)2| 



(10) 



One should note that the formula obtained here applies to 
all soliton-like objects consisting of a field configuration 
with some proper modifications. 



Now that we have developed a method to identify the 
location and velocity of strings at each time, the La- 
grangian evolution of the string network can be traced 
just as in the case of numerical simulation of the evo- 
lution of local strings based on the Nambu-Goto action. 
In fact, our simulation, being based on scalar field equa- 
tions, contains even more information, that is, we can 
find fate of intersecting string segments, namely, whether 
they reconnect or simply pass through each other, with- 
out assigning the probability of reconnection by hand 
unlike in the case of simulations based on the Nambu- 
Goto action. We can therefore calculate formation rate 
and spectrum of string loops without ambiguities. 

Thus we can calculate all the terms of eq. (|l|) directly 
from simulation data except for the last term, which can 
be evaluated by use of this equation itself. As pointed out 
in we should use ^ = 7/1^ as the energy per length 
to calculate the energy density of long strings, poo- Here 
7 is the average Lorentz factor and /i^ is the line density 
of a static string given by /is ~ 27r(T^ ln(t/(5^^/^)). 

Figure |l| depicts time evolution of ^ with three dif- 
ferent identification methods of strings. Filled squares 
represent results of our new identification scheme, while 
blank circles and blank squares correspond to our previ- 
ous method [^,^ and Vachaspati-Vilenkin algorithm, 
which was adopted in ||2^, respectively. Both methods 
estimate total string length by simply counting the num- 
ber of boxes penetrated by a string. As is easily seen, our 
previous method slightly overestimates ^ and Vachaspati- 
Vilenkin approach overestimates ^ by a factor of 1.4. We 
confirm that ^ becomes constant after some relaxation 
period and find that the more correct value of the scal- 
ing parameter in the radiation dominated era is given by 
i ~ 0.80. 

Figures ^ depict time evolution of average velocity, av- 
erage velocity squared, and average Lorentz factor. The 
first two quantities are also found to relax to constant 
values after some relaxation time, namely, {v) ~ 0.65 
and (w^) ~ 0.50 > {v) . The average Lorentz factor, 
however, has a large scatter in time, although long-time 
average is fairly constant with 7 ~ 1.8. Since string seg- 
ments moving with a speed close to light velocity have 
extremely large Lorentz factor and push up the average 
dramatically, fluctuation in the number of such string 
segments results in such large scatter. This also explains 
why the average Lorentz factor 7 = 1.6—2.0 is larger than 

1 / a/1 — (u^) and 1 / \Jl— (w)^. Thus the energy per unit 
length of a string is enhanced by a factor 7 ~ 1.8. 

Calculating the length of all infinite strings and that 
of loops at each time step and comparing with those at 
the preceding time step, the loop production parameter 
is found to be c = 0.43 — 0.53. Then taking account of 
the relation (||), the emission parameter k is calculated 
as K = 0.03 — 0.13. We therefore conclude that the di- 
rect emission of NG bosons is a subdominant channel of 
energy loss of long strings even in the case cosmic hori- 
zon scale is not much larger than the string width due to 
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the limitation of the dynamic range of numerical simula- 
tions. Thus, even if k might become significantly smaller 
in the cosmological situation because it is proportional 
to ln(i/^) as pointed out in |^,^, the scaling param- 
eter f would not increase to coincide with that of local 
strings. The result of our simulation with ln{t/S) ^ 5 
implies that k ^ 0.01 in the cosmological situations. In- 
serting this new value of k into the relation (H) together 
with c, K, and (v'^) yields ^ = 0.9 — 1.4 in the cosmologi- 
cal situations, which is still significantly smaller than the 
case of local strings. Hence, it is confirmed that although 
global strings relax to the scaling solution just as local 
strings, their number density is much smaller than the 
local ones. This is mainly because global strings inter- 
commute more often and (w^) is larger due to long-range 
attractive forces between strings. Finally, we apply our 
results into the constraint of the breaking scale fa of the 
Peccei-Quinn U(l) symmetry. Then, in accordance with 
Ref. m, /a is constrained as fa < (0.20 - 1.6) x lO^^ 
GeV for the normalized Hubble parameter h — 0.7 with 
C = 0.9-1.4 and 7~ 1.6-2.0. 

In summary, we have developed a new method to trace 
Lagrangian evolution of topological defects solving scalar 
field equation on three dimensional Eulerian lattices and 
given all the quantities characterizing cosmological evo- 
lution of the global string network, namely, energy den- 
sity, peculiar velocity, Lorentz factor, intercommutation 
rate, and emission rate of NG bosons. All these quan- 
tities can be obtained in Eulerian simulations only after 
establishing the method given in this Rapid Communica- 
tion, which enables us to extract Lagrangian information 
from simulations done in Eulerian meshes. Our method 
is directly applicable to situations in a variety of fields in 
physics. 
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FIG. 1. Time development of f is depicted. Blank circles 
represent time development of ^ for the identification method 
done in Refs. |^,^. Filled squares represent the new iden- 
tification method. Blank squares represent the identification 
method based on the Vachaspati-Vilenkin algorithm. 
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FIG. 2. Time development of average velocity, an average 
of velocity squared, and an average Lorentz factor of global 
strings is shown. 
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